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ABSTRACT 

Measuring star formation rates (SFRs) in high-z galaxies with their rest-frame ultraviolet (UV) continuum 
can be uncertain because of dust obscuration. Prior studies had used the submillimeter emission at 850 pm to 
determine the intrinsic SFRs of rest-frame UV selected galaxies, but the results suffered from the low sensitivity 
and poor resolution (~ 15")- Here, we use ultradeep Very Large Array 1.4 GHz images with ~ l"-2" reso- 
lutions to measure the intrinsic SFRs. We perform stacking analyses in the radio images centered on ~ 3500 
Lyman break galaxies (LBGs) at z ~ 4 in the Great Observatories Origins Deep Survey-North and South fields 
selected with Hubble Space Telescope/ Advanced Camera for Surveys data. The stacked radio flux is very low, 
0.08 ±0.15 /iJy, implying a mean SFR of 6 ± 1 1 M yr" 1 . This is comparable to the uncorrected mean UV 
SFRs of ~ 5 Mq yr" 1 , implying that the z ~ 4 LBGs have little dust extinction. The low SFR and dust extinc- 
tion support the previous results that z ~ 4 LBGs are in general not submillimeter galaxies. We further show 
that there is no statistically significant excess of dust-hidden star-forming components within ~ 22 kpc from 
the LBGs. 

Subject headings: dust, extinction — galaxies: evolution — galaxies: high-redshift — radio continuum: galax- 



1. INTRODUCTION 

Understanding galaxy evolution in the early Universe 
requires large samples of various kinds of high-z galax- 
ies. One important technique enabling selections of z > 
2.5 star-forming galaxies is the Lyman break technique 
(ICowie et all 1 1988b ISongaila et al.1 119901: iLillv et alJ 11991b 
iSteidel & Hamilton! [l 993b iSteidel et al.l 1 1 995b also see a re- 
view in Giavalisco 2002). The Ultra Violet (UV) spectrum of 
a high-redshift star-forming galaxy usually exhibits a Lyman 
continuum discontinuity at 912 A, which is caused by absorp- 
tion of H I in the stellar atmosphere of massive stars, and in- 
terstellar and intergalactic media. By using this feature, large 
samples of Lyman break galaxies (LBGs) can be selected in 
various redshift ranges by searching for sudden brightness 
dropouts between two adjacent broad-band images. 

An important property of LBGs is their star formation rates 
(SFRs). However, determining SFRs for z > 3 systems is still 
challenging because most of the diagnostics at low redshifts 
are unavailable when light becomes dimmer and redshifts to 
wavelengths that are harder to observe. One of the most com- 
monly used method to determine SFRs at high redshift is to 
measure the rest-frame UV continu um which shifts at z > 3 
to the optical and near-infrared (e.g., Kennicutt 1998). Unfor- 
tunately, the UV continuum can be easily atte nuated by dust, 
which results in underestimation of the SFRs. iBouwens et alJ 
(2009) attempted to determine the correction factor by us- 
ing the observed correlation between the ratio of far-infrared 
(FIR) to UV flux and the UV spectral slope (Meu rer et alJ 
1999). They found that UV continuum of z ~ 4 LBGs in the 
Great Observatories Origins Deep Survey-North (GOODS-N) 
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and South (GOODS-S) fields underestimates SFRs by factors 
of ~ 3-6. However, this result is debatable since the corre- 
lation exhibit s significant scatter i n different populati ons of 
galaxies (e.g. JCortese et al. 112006b iHowell et al. II2010I) . 

A more direct way to estimate the SFRs is to deduce them 
at longer wavelengths where light is not affected by dust ob- 
scuration. Efforts had been made using the submillimeter 
850 /im emission to determine the intrinsic SFRs in LBG s 
dPeacock et alJ2000blChapman et alJ2000HWebb et al.l2003l) . 
However, these results have low signal-to-noise ratios (S/Ns), 
presumably because the submillimeter single-dish maps were 
not deep enough, a consequence of instrumental, sky, and con- 
fusion noises (i.e., uncertainties contributed by nearby bright 
sources or faint undetected sources). In addition, an assump- 
tion of dust temperature is required to estimate the total in- 
frared (IR) luminosity. These all make the estimate of SFRs 
based on submillimeter measurements quite uncertain. 

The radio wavelength is an alternative probe of SFRs. At 
1 .4 GHz, the radio continuum is dominated by synchrotron 
radiation from relativistic electrons produced by supernovae. 
By converting the 1.4 GHz luminosity to tot al IR luminosit y 
with the well-known radio — FIR correl ation (ICondonl 1992). 
it is possible to estimate intrinsic SFRs (Kennicutt 1998]). The 
radio — FIR correlation is fairly insensitive to dust tempera- 
ture. The angular resolution in the radio can be very high 
(~ I"- 2"), meaning very little confusion noise. Because 
of these great advantages, this radio-ba sed method has been 
used in various high-z studies (e.g . . iReddv & Steidell 12004b 
Wang lFalEooaiCarilli et aTJl2008b lPannella et al] |2009h . 

The Very Large Array (VLA) provides high resolution in 
the radio, but its sensitivity is insufficient for directly detect- 
ing normal star-forming galaxies at z > 1.5. Therefore, stack- 
ing radio fluxes of large samples of LBGs to determine their 
mean SFRs is a necessary approach. IReddv & Steidell (120041) 
stacked z ~ 2 LBGs and f ound the dus t correc tion to be ~ 4.5. 
Using a similar method, Car illi et al.l ([2008) determined the 
dust correction for z ~ 3 LBGs to be 1.8. 

In order to determine the intrinsic SFRs for z ~ 4 LBGs, 
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we performed radio stacking analyses in the GOODS- 
N and GOODS-S fields. The GOODS Hubble Space 
Telescope/Advance Ca mera for Surveys (HST/ACS) data 
(IGiavalisco et al.ll2004h were used to locate ~ 3500 LBGs in 
the fields, and deep V LA 1.4 GHz images (IMiller et alj|2008t 
iMorrison et al.l l2010) were used to measure and stack the ra- 
dio fluxes. We note that although radio stacking of z ~ 4 
LBGs had also been attempted by ICarilli et ail (120081) . we 
used much deeper optical images (~ 2 mag deeper at V band 
) and radio images (~ 15% -40% deeper). We are therefore 
probing much cl oser to the typic al members in the z ~ 4 LBG 
population, as in lBouwens et al.l ((2009). 

In this paper, we first describe the observational data in 
Section [2] In Section [3] we describe our methods of select- 
ing LBGs, stacking analysis, and estimations of SFRs using 
radio and UV. Finally, we compare the two different SFRs 
and discuss the implications in Section HI and give a sum- 
mary in Section [5] We assume fio = 0.3, fl\ = 0.7, and 
Hq = 70 km s" 1 Mpc" 1 . All magnitudes are in the AB mag- 
nitude system. 

2. DATA 

2.1. HST/ACS 

The HST/ACS mult iband imaging in the GOODS fields 
(IGiava lisco et al. 2004) consists of four passbands: F435W, 
F606W, F775W, and F850LP, which are referred to as £435, 
V606, ! 775' an d ^850, respectively. We adopt the v2.0 source 
catalogs for selecting our LB G samples. Quantities mea- 
sured with SExtractor dBertin & Arnoutsl 119961) "automatic 
aperture" (MAG_AUT0, FLUX_AUTO, etc.) are used to ap- 
proximate total values. 

2.2. VLA 

The VLA 1.4 GHz images of the GO ODS-N and GOODS- 
S fie lds are descr i bed in more details in Morr ison et al.l (120101) 
and IMiller et all (120081) . respectively. In brief, the North 
field has an rms noise of ~ 4 /iJy beam" 1 at the center, and 

< 6 /iJy beam" 1 at the edges of the ACS fields. The South 
field has an rms noise of ~ 6 /iJy beam" 1 at the center, and 

< 7.5 fiiy beam" 1 at the edges of the fields. These sensitivities 
have been corrected for the primary beam response. The beam 
FWHM of the North field is 1 "7 x 1 ."6, which corresponds to 
1 1 .8 kpc x 1 1 . 1 kpc at z = 4. The beam FWHM of the South 
field is 2."8 x 1"6, which corresponds to 19.5 kpc x 11.1 kpc 
at z = 4. Data Release 2 of the South field is used here. 

3. METHOD 

3.1. Sample Selection 

"B-dropout" galaxies at z ~ 4 are selected utilizing the red- 
shifted Lyman break located between B^j and Vaae- We 
adopt the well establi shed criteria (e.g., Bec kwith et alj|2006t 
iBouwens et al.ll2007l) : 

B435 - V606 > 1 ■ 1 , 

B435 - V606 > ( V606 - Z85()) +1-1, 
V606-Z850 < 1.6, 

S/N(V m ) > 5, and S/N(i 115 ) > 3. 

In addition, compact objects (SExtractor stellarity indices 
greater than 0.8) with 1775 < 26.5 are rejected from our sam- 
ple to prevent stellar contamination. A color-color diagram 



illustrating the sample selection is shown in FigureQ] In total, 
we selected 1778 and 1679 B-dropouts in the North and South 
fields, respectively. These selection criteria efficiently select 
LBGs between z ~ 3 and z ~ 4.5, with a m ean redshift of 
3.8 (see more details in lBouwens et al.ir2007l) . The number of 
LBGs selected by us is very similar to that in lBouwens et al] 
(120071) . who used similar GOODS ACS data. 
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FIG. 1 . — A color-color diagram showing the B-dropouts selection criteria. 
Objects in the GOODS fields are plotted wi th do ts. Intersection of the three 
color selection criteria described in Section [3TT1 are shaded in grey. The se- 
lected LBGs are shown with red dots. Note that some of the objects inside 
the grey region are not selected as LBGs because they do not meet the other 
criteria (e.g., low S/N or too compact). 

3.2. Stacking Analyses 

We measured the radio fluxes of the z ~ 4 LBGs with aper- 
ture photometry at their optical positions. Small offsets (< 
0".3) between the coordinates of the radio and optical images 
had been corrected utilizing their source catalogs. Apertures 
of two different sizes were used here: small-apertures with 
radii equal to the beam FWHMs (~ 1 1 kpc at z = 4); large- 
apertures with radii equal to twice the beam FWHMs (~ 22 
kpc at z = 4). Here we first focus on the small aperture re sults 
and we will discuss the large aperture results in Section 14.41 
Elliptical apertures were used for the South field to match its 
beam shape. Fluxes measured with our aperture photome- 
try were calibrated w ith the radio catalogs ( Mille r et alJl2008b 
Morri son et al.l l2010). For sources in the radio catalogs, we 
measured their fluxes with our aperture photometry and com- 
puted the median of the ratios between our fluxes and the cat- 
alog fluxes. The ratios, which were found to be very close to 
unity for the aperture sizes we used, were applied back to our 
radio fluxes. 

Sources with radio fluxes > 100 pJy were rejected from 
our measurements because their radio fluxes do not likely re- 
flect their true SFRs. There are eight such sources. If they 
are normal star-forming galaxies at z = 4, their IR luminosity 
would be > 5 x 10° L (or SFR > 7000 M yr" 1 ), which 
is unusually large. We inspected these sources individually. 
They either are affected by nearby bright radio sources that 
do not appear to be at the same redshifts of the LBGs, or do 
not appear to be associated wi th known submill i meter galax- 
ies (SMGs) in the samples p f lWang et al.l d2004l).lPerera et alj 
(2008), iDevlin et al.l d2009l) . and lWeifiet alj ( 120091) . There- 
fore, they are highly unlikely high-redshift ultraluminous star- 
bursting galaxies. To avoid the bias from unrelated nearby 
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bright sources and from active galactic nuclei (AGNs), we did 
not include these eight sources in our stacking analyses. For 
those sources with radio fluxes < 100 /iJy, we then averaged 
their radio fluxes and subtracted a background value (see be- 
low) from the means to get the final stacked radio fluxes. 

The uncertainties of our stacked fluxes were estimated with 
Monte Carlo simulations, with an assumption that the galax- 
ies are distributed randomly over the map. We measured the 
mean radio fluxes at random positions, and the number of ran- 
dom positions is the same as that of the < 100 fjJy B-dropouts. 
The same apertures were used and the same rejection crite- 
rion of < 100 /iJy was adopted. This measurement is re- 
peated 10,000 times, and the mean radio fluxes have a fairly 
Gaussian-like distribution. We then calculated the mean and 
the dispersion of these 10,000 measurements. The mean was 
subtracted from the stacked radio flux of the LBG sample to 
form the final stacked radio flux. This subtraction is to ac- 
count for the effect of imperfect clean and chance projection 
of random radio sources in our flux apertures. Likewise, the 
dispersion can then represent the uncertainty of the stacked 
radio flux. 

There are two caveats in the above procedures. First, if 
Z ~ 4 LBGs are clustered at scales similar to the sizes of 
our apertures, our assumption that LBGs are randomly dis- 
tributed would break down. Our stacking method would then 
overestimate the radio flux es. We investig a ted clu stering with 
the same method used by Marsde n~et al.l (12009), We found 
that the variance-to-mean ratio of source numbers inside ran- 
domly placed r < 10" circles of our sample exceeds that of 
random distribution by less than 0.1, and becomes smaller 
with smaller radii. Therefore, we conclude that there is no 
significant small-scale clustering in our LBG sample and the 
flux overestimate caused by small-scale clustering can be ne- 
glected. 

Second, the faulty estimate of uncertainties can be caused 
by concentration of sources in high or low sensitivity regions 
(i.e., the slightly uneven sensitivity distribution caused by the 
primary beam falloff in the GOODS-N or by the mosaicking 
in the GOODS-S). Instead of stacking at random positions to 
estimate the uncertainty, we also carried out the simulation 
based on the source positions of the real LBG sample. We 
added an offset to the positions of the LBGs and stacked their 
fluxes. We repeated this for 10,000 times with different ran- 
dom offsets of 15"— 60", and calculated the mean and the dis- 
persion. The results are within 30% to those measured from 
random positions, suggesting that LBGs are not concentrated 
in specific regions. 

With the above methods, we measured the final stacked 
fluxes and flux errors of the LBGs in the two fields separately. 
For each aperture size, we combined the two final stacked 
fluxes by weighting them with the inverse of the square-errors 
to form the combined flux. The results are summarized in Ta- 
ble Q] The stacked radio signal is very weak, consistent with 
zero within the noise. 

3.3. Star Formation Rate 
3.3.1. Radio SFR 

The stacked radio fluxes can be converted to SFRs. As- 
suming a universal synchrotron emission spectral index of 
a = -0.8, we convert the stacked radio fluxes S sta ck to the rest- 
frame 1.4 GHz luminosity densities Li.4ghz, i-e-, 



TABLE 1 
Mean Radio Fluxes 



Field 


Small Aperture" 


Large Aperture b 




^LBG 


^random 


^LBG 


^random 




(MJy) 


(MJy) 


OuJy) 


(MJy) 


GOODS-N 


-0.05 ±0.1 8 


0.36 ±0.1 8 


0.29 ±0.28 


1.10 ±0.28 


GOODS-S 


0.36 ±0.26 


0.10 ±0.26 


0.54 ±0.54 


-0.03 ±0.54 


Combined 


0.08 ±0.15 




0.34 ±0.25 





NOTE. — ^random is the mean radio flux of random positions derived with 
the Monte Carlo simulation in Section [3~2l This mean is subtracted from the 
measured flux of the LBG sample, and its error is propagated to the mean of 
the LBG sample. Flbg listed has already been subtracted by F mn( i om . 



"Apertures with radii equal to the beam FWHM. 
b Apertures with radii equal to twice the beam FWHM. 



where d\ is the luminosity distance and z is the m ean redshift 
of 3.8. With the loc al radio — FIR correlation dHelou et alj 
119851; ICondonl [19921) . L lAGVa can be converted to FIR lu- 
minosity Lfir(40-i20), approximately the total luminosity be- 
tween A = 42.5 /im and A = 122.5 fim, with the following re- 
lation: 



q = log 



^FIR(40-120) 

3.75 x 10 12 W 



-log 



M.4GHz 

WHz" 1 



(2) 



with q = 2.34 ± O.Ol dYun et al. 1120011) . Although this value 
was derived locally (z < 0.15), recent s t udies suggest weak 
evoluti on out to z ~ 2 (Sar gent et al. 1 120 10t llvison et aTl 
l2010alfbT) . Extrapolating q to z = 4 using the relat ion q oc 
(1 +zV with 7 = -0.04 ± 0.03 (llvison et al. 1 1201 Obi) implies 
that Lfir (and therefore radio SFRs) will be lowered by only 
~ 30% ± 20%, which is insignificant compared to the uncer- 
tainties in our results. 

SFRs can th en be estimated with the conversion in 
!Kennicuttl(fT998h : 



SFR 



CM© yr" 1 ) 



= 4.5 x 10 



-44MR(8-1()00) 



erg s" 



(3) 



£l.4GHz = 47rt/ / 2 5 stac k(l +Z)~ 



a+a) 



(1) 



where LiR(g_iooo) is approximately the total luminosity be- 
tween 8 /.tm and 1000 /im. This estimation is within ~ 30% 
to other published calibrations. We derived the ratio between 
£fir(40-i20) and Lir(8-iooo) by computing these two quantities 
on model spectral energ y distributions o f six nearby normal 
and starburst galaxies (ISilva et all 119981) . using their origi- 
nal definitions over the Infra red Astronomical Satellite (IRAS) 
bands (see the summary in Sanders & Mirabel 1996). The 
£ir(8-iooo) to L F ir ( 4o-i20) ratio has a range of 1.71-2.32, with 
a weak anti-correlation between the ratio and IR luminosity 
and a mean of 2.05. This m ean ratio is slightly higher than 
that used in lYun et al. I (12001b of 1.5, which is based on mea- 
sure ments in luminous IRAS galaxies and s t arburst galax- 
ies dSandersetal.il 199 It iMeurer et all 119991; iCalzetti et all 
2000), but is closer to the values in nearby normal spiral 
galaxies and low luminosity starbursts (~ 2 for M 51 and 
M100, and 2.3 for M 82, based on the Silva et al. templates). 
Using this value, we derived the SFRs listed in Table|2j which 
are consistent with normal galaxies or low-luminosity star- 
bursts. 

3.3.2. UVSFR 

The (775 and Vw6 bands correspond to rest-frame UV, and 
can be used to calculate SFRs, uncorrected for extinction 
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TABLE 2 

Star Formation Rates Derived with Various Fluxes 



Radio Optical 
Small aperture Large aperture (775 Z850 
SFRs(M yr') 6.0 ±11.0 25.7 ±18.8 4.92 5.14 



(lMadauetal.ll 19981) . with 

r Luv(ergs s" 1 Hz -1 ) 
SFR(M yr" 1 ) = S g , (4) 

where Luv is the UV luminosity density at the same redshift 
of3.8,andCis8.0x 10 27 and 7.9 x 10 27 at 1500Aand2800A, 
respectively. We list in Table [2] the mean UV SFRs of the B- 
dropouts with radio fluxes < 100 fiJy. The uncorrected SFRs 
derived from / 775 and Va)b are remarkably similar to each other 
as well as to that derived from the radio stacking analyses, 
indicating very little extinction in the observed rest-frame UV 
emission. 

4. RESULT AND DISCUSSION 

4.1. Dust extinction in z ~ 4 LBGs 

In Figure [2] we plot the z&so magnitudes versus the small- 
aperture radio fluxes and the corresponding SFRs of all the 
Z?-dropouts in the two fields, with the red points indicating 
the mean (< 100 /iJy sources). As can be seen in the fig- 
ure and Table the average SFRs deduced from the radio 
and the rest-frame UV are consistent within measurement un- 
certainties. The relatively small SFR derived from the radio 
flux suggests little extinction in the z ~ 4 LBGs. Unfortu- 
nately, although this is by far the de epest radio stackin g on 
Z ~ 4 LBGs (~ 3 X deeper than that in lCarilli et alll2008l) . the 
error in the stacked radio flux is still too large to pin down 
the extinction correction. An extinction correction of > 6 is 
ruled out at 2. 2 a and an extinction correction of > 4 is ruled 
out at 1 .3 a. iBouwens et al.l (120091) derived extinction cor- 
rections of ~ 3 - 6 for z ~ 4 LBGs based on their UV con- 
tinuum slopes. Our result slightly favors a lo wer value and is 
margin ally (~ 2 a) consistent with the result of Bo uwens et alJ 
(2009). 

We note that an even stricter upper limit on the radio 
SFR can be placed if any of the following is true: (1) the 
radio — FIR corre latio n q factor exhibits a decline with red- 
shift (see Section l43T >: or (2) the conversi on from Lfir(40-120) 
to Lir(8-iooo) is overestimated (see Section l3.3.U : o r (3) there 
is AGN contribution in our radio SFR (see Section FOI i. The 
combined effect of all the above may not be negligible and 
can make the extinction correction for z ~ 4 LBGs even lower. 
H owever, this remain s to be tested by future observations. 

ICarilli et"afl ([2008) found a low extinction correction of 1 .8 
on z ~ 3 LBGs. On the other hand, on z ~ 4 LBGs, they had 
a 2 a detection of 0.83 ±0.42 /iJy on 1447 Z?-dropouts, which 
also suggests a low extinction correction. This radio flux is 
much higher than ours (0.08 ±0.15 /iJy) but their ground- 
based V-band limiting magnitude is also much shallower than 
ours (by ~ 2 mag). In other words, we are probing different 
LBG lu minosities. Compar ing our z ~ 4 result with the z ~ 4 
result in lCarilli et al.l (120081) is therefore not straightforward. 

4.2. AGN contamination 

To avoid AGN contamination, a flux cut of < 100 /iJy was 
set and eight sources were rejected (one in the north and seven 
in the south). If these sources comprise all the AGNs in our 
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FIG. 2.— HST zsso magnitude vs. 1.4 GHz VIA flux of the ~ 3500 B- 
dropouts in the GOODS fields. Orange symbols indicate objects not included 
in our stacking analyses. The radio fluxes are measured with small apertures 
with radii equal to the beam FWHMs (~ 1 1 kpc at z, = 4). The z,gsa uncertainty 
of each data point is indicated with vertical bars. The radio uncertainties, 
which depend on the locations on the maps, are not shown to avoid confusion. 
Instead, we plot the typical radio flux uncertainties for galaxies close to the 
center of the 2 fields at the lower-right corner for reference (upper: North, 
lower: South). The inferred SFR corresponding to the measured flux density 
at 1.4 GHz (see text for details) is shown on the upper x-axis. The average 
of all the black data points are shown with the red symbols, with the radio 
uncertainties estimated from the MonteCarlo simulation. The black curves 
indicate the locus of equal SFR deduced from the z%ia and the radio. The 
color coded curves, from left to right, indicate 1/6, 1/4, 1/2, 2, 4, and 6 times 
dust corrections from the UV to radio SFR. A blow up to the stacked data 
point is shown as an insert. 

sample, our sample would have a very low AGN fraction of 
~ 0.2%. For comparison, the AGN fraction in LBGs at z ~ 
3 was estimated to b e ~ 3%, based on optical spectroscopy 
( Shapl ev et al. 1 120031) . This difference could be due to the 
fact that our LBGs are much less luminous than those in the 
spectroscopic sample. 

It is possible to roughly separate AGNs and star forming 
galaxies based on the radio power. In the lo cal radio lumi- 
nosity function (e.g.. iMauch & Sadler 1 120071) . the 1.4 GHz 
power that divides AGNs and starbursts is ~ 10 23 W Hz" 1 . 
ICowie etail (120041) show that this AGN/starburst division in- 
creases to ~ 10 24 W Hz" 1 at z ~ 1. This power corresponds 
to ~ 10 t iSy for z = 4. How ever, the increasing trend ob- 
served by Cowie et al. (2004) can continue to z 1, moving 
the AGN/starburst division to a radio flux greater than 10 /iJy. 
Furthermore, high-redshift flat-spectrum radio sources (pre- 
sumably AGNs) have a more negative ^-correction than that 
of steep-spectrum sources (starbursts). This also increases 
the apparent flux division between AGNs and starbursts at 
observed-frame 1.4 GHz (by roughly a factor of 2). We thus 
expect sources below our 100 fjJy cut to be dominated by star 
forming galaxies. 

It is interesting to test what happens if we assume that > 50 
/iJy sources are also AGNs and remove them from the stack- 
ing analyses. To do this, we rejected 13 sources that are 
brighter than 50 /iJy and introduced the same 50 /iJy cut 
in the Monte Carlo simulation. We found a stacked flux of 
0.05 ± 0. 14 /iJy, corresponding to a mean SFR of 3.8 ± 10.5 
M Q yr" 1 . The stacked flux decreases further if we keep lower- 
ing the flux cut but the change is well within the noise. Never- 
theless, it is obvious that we will get a lower mean SFR if we 
assume more AGN contribution. Although we are limited by 
noise here, our conclusions that z ~ 4 LBGs have low SFRs 
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and the mean extinction correction seems to be lower than 5 
are thus still valid. 

4.3. The Radio — FIR Correlation 



iCarilli etail ((2008) discussed the use of the local radio — 
FIR correlation to derive SFRs of high-redshift galaxies. They 
mentioned several possibilities that this correlation may be 
different at high redshift. I ndeed, a weak evolution has been 
sugges ted by recent studies ( Sargen t et al. l2010lHvison et al. I 
2010a,b). However, as mentioned in Section [3.3. II even if we 
extrapo late the evolution based on a declining q (llvison et al. I 
2010b) and adopt a Lir(8-iooo) to Lfir(40-120) ratio that is 25% 
lower, our main conclusions would not change and this would 
make the SFR of the LBGs even lower. Furthermore, the local 
radio — FIR correlation is only used to calibrate the relation 
between SFR and radio power. Even if there is an evolution 
in the radio — FIR correlation at high redshift, we cannot be 
certain that the evolution arises from the radio part of the cor- 
relation and affects the radio — SFR relation. In short, we see 
neither evidence nor strong argument that the radio SFR is 
significantly biased by the assumption of the local radio-FIR 
correlation. 

4.4. Implications for SMGs 

The relatively low dust correction and SFRs imply that 
Z ~ 4 LBGs in general do not have extreme star formation 
activities enshrouded in dust, unlike that in SMGs. If the to- 
tal SFR of our z ~ 4 LBGs is only contributed by SMGs with 
SFR of 1000 M yr _1 and if other non-SMG LBGs have no 
star forming activities, then there can be no more than 21 ± 39 
SMGs in our L BG sample. In comparison, the source counts 
at 850 fim (ICoppin et alj2 006) suggest ~ 350 SMGs in a field 
of the combined GOODS size, with 5^50 > 1.4 mJy and a 
mean SFR of 1000 M Q yr" 1 . This means that even if the total 
SFR of our LBG sample is dominated by z ~ 4 SMGs, they 
can only account for < 10% of the entire SMG population at 
$850 > 1 .4 mJy. This implies that either LBGs and SMGs have 
little overlap at z ~ 4 or most SMGs are at z -C 4. We note 
that previous studies in th e submillimeter showed that LBGs 
are in general not SMGs ((Chapman et al. 20061 iPeacock et all 
120051 IWebb et al.ll2003h. an d that most SMGs are at z < 3.5 
dChapman et al.ll2003L 12005). Our results are consistent with 
these. 

Another scenario related to SMGs worth investigating is 
dust-hidden companions sitting cl ose to the LBGs. A n ex- 
ample of this is the z = 4.5 SMG in lCapak et ail (120081) . This 
galaxy is independently selected as an LBG, and has an ul- 
traluminous dusty component showing up at > 3.6/im next to 
its rest-frame UV emission. To know whether similar systems 



are common at z ~ 4, we performed the same radio stacking 
using the large apertures with radii twice the beam FWHM 
(~ 22 kpc at z = 4) to measure the radio fluxes. As shown in 
Figure[3]and Table[2] the radio SFR is slightly higher but still 
consistent with the UV SFR, implying that there is no statis- 
tically significant excess of dust-hidden star-forming compo- 
nents close to LBGs. 

5. SUMMARY 

We employed radio stacking analyses on ~ 3500 LBGs at 
z ~ 4 in the GOODS-N and GOODS-S fields. The mean ra- 
dio flux was converted to FIR luminosity via the well-known 
radio — FIR correlation, and then converted to an intrinsic 
SFR, which is 6 ± 1 IMq yr" 1 . By comparing with SFRs esti- 
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FIG . 3 . — Same as Figureff] but the radio fluxes are measured with apertures 
of radii twice the beam FWHM (~ 22 kpc at z = 4). 

mated from rest-frame UV, we found a maximum UV extinc- 
tion correction of ~ 6. This is roughl y consistent with tha t 
derived from UV continuum slopes bv lBouwens et al.l d2009l) 
with a very similar sample. The low extinction correction con- 
firms that z ~ 4 LBGs are in general not SMGs. We also in- 
vestigated the possibility of finding dust-hidden companions 
close to the LBGs at z ~ 4. The stacked radio flux within 
~ 22 kpc to the LBGs implies a mean SFR that is still con- 
sistent with the UV inferred SFR, which suggests no statis- 
tically significant excess of dust-hidden star-forming compo- 
nents close to LBGs. 
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